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EFFECT  OF  MAGNETIC  SHEAR  ON  THE  ELECTROSTATIC  ION  CYCLOTRON 
WAVES  IN  MAGNETOSPHERIC  PLASMAS 

I.  INTRODUCTION 

The  current  driven  ion-cyclotron  instability  (CDICI)  has  been  of 
considerable  interest  to  space  plasma  physicists.  Recent  observations  by 
Kintner  (1981)  and  Kintner  et  al.  (1978)  have  made  it  very  topical.  Drummond 
and  Rosenbluth  (1962)  were  first  to  examine  this  instability  analytically 

while  Kindel  and  Kennel  (1971)  studied  it  in  the  context  of  space  plasmas  in 
the  earth's  magnetosphere.  The  results  from  Kindel  and  Kennel  have  frequently 
been  applied  to  observations  involving  the  electrostatic  CDICI  in  space.  It 
should  be  noted,  however,  that  the  analysis  of  both  Drummond  and  Rosenbluth 
and  Kindel  and  Kennel  is  a  local  analysis.  They  consider  a  uniform  zeroth 
order  magnetic  field  (i.e.,  B_  *  BQ2),  neglecting  the  self-consistent  magnetic 
field  generated  by  field  aligned  currents.  This  magnetic  field  (usually 

small)  will  give  rise  to  a  shear  in  the  zero  order  magnetic  field  (see  Figure 
1)  and  consequently  make  the  zero  order  field  space  dependent  (i.e.,  JS(x)  = 
B^x)y  +  Bzz,  where  shear  is  in  the  x-direction) . 

Magnetic  shear  in  general  is  a  damping  agent  and  can  significantly  alter 
the  local  mode  structure  (e.g.,  it  was  shown  by  Ross  and  Mahajan,  (1978)  that 

an  infinitesimal  shear  can  completely  damp  the  universal  drift  instability). 

Recently  Waltz  and  Dominguez  (1981),  motivated  by  the  TFR  (1978)  experiment, 
have  provided  numerical  results  for  the  behavior  of  CDICI  in  a  sheared 
magnetic  field  pertinent  to  the  TFR  (1978)  parameters.  Ganguli  and  liakshi 
(1981)  have  given  a  detailed  analytical  treatment  of  the  CDICI  in  a  sheared 
magnetic  field  and  have  concluded  that  even  a  small  shear  can  give  a 
significant  reduction  in  the  growth  rates.  It  also  greatly  reduces  the  region 
of  unstable  perpendicular  wavelengths. 

Manuscript  submitted  March  29,1982. 


1 


In  this  brief  report  we  illustrate  the  importance  of  the  non-local 
treatment  of  the  CDIC1  for  the  space  plasmas.  We  shall  show  that  the  nonlocal 
theory  leads  to  enhanced  coherence  in  the  electrostatic  CDICI,  consistent  with 
the  observations  of  Kintner  (1981). 


II.  THEORY 


a.  Assumptions 

We  consider  a  slab  geometry  where  both  the  electrons  and  the  ions 
have  a  Gaussian  distribution  function.  The  plasma  has  negligible 
MM  nkT/8*  B2  «  1)  and  the  electrons  drift  with  respect  to  the  ions  along 
the  magnetic  field  with  a  velocity  V^.  both  species  have  a  finite  non-zero 
temperature  T^  =  m^v^2/2  kfi,  or  denoting  the  species,  is  the  mass,  is 

the  thermal  velocity  and  kg  is  the  Boltzmann's  constant.  We  neglect 

(kX^)2  <<  1  where  k  is  the  wavevector  in  arbitrary  direction  while  is  the 

Debye  length.  Since  for  most  space  plasmas  of  interest  T^  ^  T.  we  shall 
restrict  the  analysis  here  for  the  temperature  ratio  x  *  T^/T  ^  1.  The 
ambient  magnetic  field  is  in  the  £-2  plane  and  contains  a  shear  in  the  R- 

direction.  This  field  is  given  by 

B  *=  Bq{z  +  (x/Lg)y}  ,  (1) 

_1  _1 

with  x/Lg  <<  1  and  Lg  =  (99/3x)  where  0  =  tan  (By/Bz).  Thus  Lg  is  the 

scale  length  characterizing  the  magnetic  shear.  In  the  absence  of  shear  the 

field  configuration  is  »  BQ2.  P^asma  Just  described  is  unstable  to  the 

CDICI  when  u>  -  k^V^  <  0  for  w  and  kfl  consistent  with  the  dispersion 

relation.  The  waves  at  maximum  growth  (for  T  *  T.)  are  characterized 

e  i 

by  (o^  ~  ft^,  ft^  =  eB  /m^c,  the  ion  Larmour  frequency;  Y  «  ;  and  kp^  ®  0(1) 

°  1/2 

where  p.  *=  (2k^T./m,)  /ft.,  is  the  ion  larmour  radius.  Electron  Landau  and 

l  b  l  i  1 

ion-cyclotron  dampings  reduce  the  growth  rate  depending  on  the  magnitude 
of  (see  Ganguli  and  Bakshi  (1981)). 

The  above  description  of  the  plasma  gets  significantly  modified  by  the 
introduction  of  a  shear  in  the  magnetic  field  as  given  in  (1).  The  magnetic 
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field  rotates  in  the  £-2  plane  (see  Figure  1)  as  a  function  of  x.  If  at  x  =  0 
we  have  *  0  then  at  x  -  XjWe  see  that  k^  *  0  since  at  xj  the  magnetic  field 
has  rotated  by  an  angle  0*  =  xi/Lg*  Thus  the  dispersive  properties  of  the 
plasma  are  also  a  function  of  x.  It  is  this  fact  that  changes  the  boundary 
condition  (from  plane  waves  to  outgoing  energy  flux  condition)  which  the  local 
theory  fails  to  account  for,  thereby  giving  dubious  results.  We  introduce 
shear  (i)  locally,  by  replacing  by  k°  +  sk^x,  where  s  =  1/Lg  and  (ii) 

globally  by  replacing  ikx  by  M  M#  Since  the  magnitudes  of  the  magnetic 
shear  of  interest  to  us  is  quite  small,  we  neglect  the  orbital  effects  of 
shear  [Bakshi,  Bellew,  Ganguli  and  Satyanarayana  (1977),  Bellew  (1978), 
Ganguli  (1980)  and  Linsker  (1981)]  arising  mainly  out  of  the  shear  kinematic 
drift  of  the  particles  in  the  sheared  magnetic  field* 
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b.  Dispersion  Relation 


The  general  dispersion  relation  for  CDICI  in  the  absence  of  shear  is 

given  by 


SI  [i  +  (— 

1  |k«,Va 


u-nfl  -k„V, 

r _ 2 _ Li°n,0 

1  »VV«  JJ  ’ 


(2) 


where  =  ^2p2/2,  ^(b)  =  In(b)e  ^  and  In(b)  are  modified  Bessel's  function. 

As  previously  described  we  study  the  effects  of  the  magnetic  shear  by 
replacing  kfi  by  kyu,  where  u  *  sx(  =  k(|(x)/k^)  and  ikx  by  —  in  equation  (2) 
(we  have  assumed  *  0).  Transforming  to  the  ion  frame  (i.e.,  setting  Vdi  * 
0),  retaining  the  n  =  0  term  for  the  electrons  and  n  *  0,  +  1,  +  2  terms  for 
the  ions  (since  we  wish  to  study  the  first  harmonic  only)  and  expanding 
rn(b-(pjr/2)32/ax2)  to  0(32/3x2)  (here  b  -  k2p2/2)  we  obtain  a  second  order 
differential  equation 

{-A(x)  Ijjr  +  Qi(x)}$(x)  =  0  .  (3) 

The  equation  can  be  reduced  to  (see  Ganguli  and  Bakshi  (1981)), 

{(PjS)2  |^7  +  Q(u)}$(u)  ■  0,  (4) 


where. 


Q(u) 


Ql(u) 

~aoT)  » 


Qi (u)  -  l  +  t  +  ciroz<c1)  +  c1r1{z((i-P)ci)  +  z  ((l+p)^)} 
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+  +  z((l+2p)C1)} 


k„V,  k a  Vj 

+  1  «.(*  - -r1  W'.«- -/  »• 

a(u)  -  )[z(ct )r'  +  {zta-pkj))  +  z((i+PKi)}r{ 
+  {z((l-2p)ci )  +  z((l+2p)c1)}r|], 


and 


ar 

Jb 


We  have  assumed  k£  «  k*  in  neglecting  higher  order  derivatives  in  the 
expansion  of 

Following  Ganguli  and  Bakshi  (1981)  we  expand  the  ’’potential**,  Q(u)  of 
the  differential  equation  (4)  around  £  *  (u-uq)  =  0  to  0(£2).  Here  uQ  is  the 
angle  of  propagation  for  the  maximum  growth  rate  in  the  absence  of  shear.  We 
thus  obtain, 


{(p^) 


2 


+  a  +  ~y~  n2}*(n)  *  o, 


(5) 


where. 


n  -  1  +  Q'/Q 

o  o 


‘-%-X  2/2V' 
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The  dispersion  relation  the  under  outgoing  energy  flux  boundary  condition 
similar  to  Pearlstein  and  Berk  (1969)  is, 


1  /2 


Q(io,u  )  =  (2  4+l)(p  .$)(-Q  "/  2)  '  +  Q'  ^/2Q  "  (6) 

O  1  O  0  0 


Note  Q(w,u^)  =  0  yields  the  local  dispersion  relation  maximized  over  k ^  , 
while  the  right  hand  side  of  the  equation  (6)  arises  out  of  the  non-local 
treatment.  The  first  term  in  the  right  hand  side  is  proportional  to  the 
magnetic  shear  and  vanishes  in  the  zero  shear  limit.  The  second  tern,  is  no< 
explicitly  shear  dependent.  This  term  can  contribute  significantly  ev-tu  r.ir 
i nf I nl tesimal  shear  thereby  making  the  nonlocal  dispersion  relation  solutions 
much  different  from  those  of  the  local  dispersion  relation. 
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III.  RESULTS 

The  non-local  dispersion  relation  (6)  is  solved  numerically  by  iterative 
technique.  We  simplify  the  "potential"  by  expanding  all  the  ion  Z-functions 
except  n=l  for  large  arguments,  retain  terms  to  0(l/£^)  and  treat  A(u )  as  a 
constant  evaluated  at  uQ.  [For  details  see  Ganguli  and  bakshi  (l^bl)J. 

Figure  2  shows  a  plot  of  Y /&  ^  against  b  for  various  values  ot  the  shear 

and  t  =  .5,  ]i  -  1837  and  V^/ve  *  .23.  Note  the  difference  between  the  local 

solution  and  the  zero  shear  limit  of  ttie  non-local  solution.  This  large 
reduction  in  the  growth  rate  is  due  to  the  second  term  in  equation  (0)  arisin0 
out  of  the  non-local  treatment.  Figure  2  also  shows  that  the  value  of  b  for 

which  maximum  growth  rate  occurs  remains  essentially  unaltered  from  the  local 

theory,  even  when  shear  is  increased.  In  addition  the  expected  reduction  ol 
the  maximum  growth  rate,  the  primary  effect  ot  increased  shear  is  to  narrow 
the  band  of  unstable  perpendicular  wavelengths.  This  result  is  responsible 
for  enhanced  coherence  of  the  wave  as  we  explain  in  the  following  section. 
The  uneven  spacing  of  the  p^s  =  0  ,  .003  and  .01  curves  shows  that  non¬ 
linearity  in  the  shear  parameter  (through  the  implicit  dependence 

of  u)  on  (p_.  s )  )al ready  sets  in  at  p^s  =  .01. 

In  Figure  3  and  4  we  provide  a  y/ti  versus  b  plots  tor  realistic 

space  plasma  parameters  relevant  to  the  S3-3  data.  In  Figure  3  we  consider 

_6 

a  hydrogen  plasma  (u  =  1837),  x  -  .73,  V^/ve  =  *3  and  (P^)  ®  10  .  Note 

once  again  the  big  drop  in  the  local  growtn  rate  as  well  as  the  sigui  t  icaut 

reduction  of  the  unstable  perpendicular  wavevector  band.  Figure  4 
illustrates  the  same  features  for  x  =  1  and  V^/ve  =  *  33.  The  consequence  or 
the  non-local  treatment  is  even  more  drastic  in  this  case.  The  nonlocal 
treatment  completely  damps  the  instability. 
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We  have  also  examined  an  oxygen  plasma  and  see  that  the  threshold  is  much 
lower  compared  with  the  hydrogen  plasma.  The  effect  of  the  non-local  term  on 
a  multispecies  plasma  is  an  interesting  question  which  is  now  being 
investigated  and  will  be  reported  elsewhere. 


ot 


Fig.  2  —  A  plot  of  the  growth  rate  (t/Hj)  against  b  for  various  shear  values.  Here  r  =  •  5, 
=  1837,  and  V<j/ve  ■  •  25.  The  big  difference  between  the  local  solution  and  the  (pjsl  * 
limit  of  the  non-local  solution  is  due  to  a  shear  independent  term  arising  from  the 
non-iocal  treatment.  The  curves  also  indicate  that  the  effect  of  the  shear  is  more  pro¬ 
nounced  for  larger  b. 


Fig.  3  —  A  plot  similiar  to  Figure  2.  Here  t  *  •  75,  n  =  1837  and  V,j/ve  =  •  3.  We  provide 
only  the  local  and  the  non-local  spectrum.  The  non-local  curve  has  a  very  small  shear,  (pjs) 
=  10"6.  Note  the  significant  reduction  in  the  unstable  perpendicular  wavevector  band. 


Fig.  4  —  Growth  rate  versus  b  plot  for  n  *  1837,  and  V«j/Ve  =  '  35.  The  effect  of  the 
non-local  treatment  is  very  drastic  in  this  case.  The  instability  is  totally  suppressed. 
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IV.  ENHANCED  COHERENCE  DUE  TO  MAGNETIC  SHEAR: 

Kintner  (1981)  and  Kintner  et  al.  (1978)  have  reported  observations  of 
the  electrostatic  CDICI  for  T  &  1  with  the  electron  temperatures  in  tile  range 
of  a  few  eV.  In  this  temperature  range  the  instability  is  observed  with  a 
high  degree  of  coherence.  Since  the  measurements  were  made  in  the  satellite 
frame  (a  frame  moving  with  respect  to  the  plasma)  the  frequencies  are  Doppler 
shifted  [Frederick  and  Coroniti  (1976)],  i.e., 

w'  *  oj  +  k* v  ,  (7) 

—  — s 

where  is  the  observed  frequency  in  the  satellite  trame,  k  is  the  wavevector 
and  Vg  is  the  satellite  velocity.  It  can  be  seen  from  (7)  that  it  the  second 
term  is  comparable  to  the  first  for  a  large  band  of  unstable  k  the  observed 
frequency  will  be  highly  incoherent.  In  the  case  ot  Si-3  data  the  satellite 
orbit  is  almost  perpendicular  to  the  magnetic  field  and  hence  k_  in  {7)  is 
just  k.  • 

Consider  Figure  3  where  T^  -  (4/3)T^f  V^/v^  *  .3  and  p  *  1837.  We  have 

plotted  ( Y/^i )  versus  b  =  (k^p^/2).  We  ignore  the  contribution  due  to  in  b 

since  k  <<  k  was  assumed  earlier, 
x  y 

By  definition  we  have. 


2b, 


thus 

ky  -  (2b)1/2  Oj/Vj. 


For  the  real  part  ot  the  frequency  Is  roughly  1.2  at  the  maximum 
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growth.  With  these  substitutions  in  (7)  we  see  that 


,..,1/2  1/2  v 

“  1.2  (1  +t!|)  v  8)  •  (8) 

c 

As  given  by  Kintner  et  al.  (1978)  we  consider  the  electron  temperature  to 
be  around  one  eV  which  means  that  the  electron  speed  is  60U  km/s  and 
using  vg  to  be  7  km/s  and  for  the  parameters  of  Figure  3  we  get, 

1/2 

ur  «  1.2  fl  (1  +  .7(b)  '). 

Thus,  at  b  ^  .8  where  the  growth  rate  maximizes  for  t  =  .75  the  second  term  is 
already  comparable  to  the  first  and  gets  larger  for  larger  b.  At  b  =  2  they 
are  both  of  the  same  order.  However  the  non-local  curve  in  the  Figure  3 
shows  that  for  the  same  parameters  the  wave  grows  in  a  narrower  wavelength 
band  beyond  which  the  second  terra  vanishes  thus  enhancing  coherence.  A 
comparison  between  the  local  and  the  non-local  curves  shows  that  the  band  of 
unstable  wavelength  is  reducd  by  more  than  half.  For  higher  shear  (see  Figure 
2)  this  band  can  be  reduced  to  less  than  a  quarter  of  the  local  value, 
thereby  giving  an  even  higher  degree  of  coherence. 
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V.  CONCLUSION 


The  effects  of  magnetic  shear  on  electrostatic  CDICI  for  the  parameters 
relevant  to  the  S3-3  satellite  data  have  been  examined.  Shear  produces  a 
large  reduction  in  the  local  growth  rate  and  significantly  reduces  the  band  of 
unstable  perpendicular  wavelengths.  The  latter  effect  leads  to  enhanced 
coherence  as  observed  in  the  S3-3  data. 

Although  we  have  shown  that  a  shear  in  the  magnetic  field  can  lead  to 
enhanced  coherence  in  the  electrostatic  CDICI,  there  can  be  other  mechanisms 
contributing  to  this  phenomenon  along  with  the  magnetic  shear.  Recently  Okuda 
and  Ashour-Abdalla  (1981)  have  shown  that  non-linear  saturation  can  also  lead 
to  coherence.  However,  it  should  be  noted  that  in  their  treatment  they  have 
used  a  plane  wave  boundary  condition  corresponding  to  the  local  analysis. 
From  our  studies  we  find  that  the  marginal  stability  can  be  signif icant ly 
altered  by  non-local  effects.  This  may  lead  to  alterations  in  the  inferences 
of  Okuda  and  Ashour-Abdalla  (1981)  and  in  other  numerical  simulations  of  the 
electrostatic  CDICI  using  a  plane  wave  boundary  condition. 

We  have  also  examined  the  effect  of  the  magnetic  shear  on  heavier  species 
and  have  found  it  to  be  equally  important.  An  interesting  situation  of 
considerable  practical  importance  is  the  effect  of  the  magnetic  shear  on  a 
multispecies  plasma.  Futhermore,  the  effect  of  magnetic  shear  on  the  Current 
Driven  Ion  Accoustic  Instability,  and  hence,  the  effect  of  shear  on  the 
relationship  of  the  critical  velocity  on  the  temperature  ratio  [Kindel  and 
Kennel  (1971)]  is  an  important  question.  These  problems  are  now  being 
investigated  and  will  be  reported  elsewhere. 
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